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, We present an alternative method for constructing the exact and approximate 

Q ' sohitions of electromagnetic wave equations whose source terms are arbitrary 

order multipoles on a curved spacetime. The developed method is based on 
^J^' the higher-order Green's functions for wave equations which are defined as dis- 

tributions that satisfy wave equations with the corresponding order covariant 
derivatives of the Dirac delta function as the source terms. The constructed 
^ I solution is applied to the study of various geometric effects on the generation 

d ' and propagation of electromagnetic wave tails to first order in the Riemann 

tensor. Generally the received radiation tail occurs after a time delay which 
represents geometrical backscattering by the central gravitational source. It is 
shown that for an arbitrary weak gravitational field it is valid that the truly 
nonlocal wave-propagation correction (the tail term) has a universal form 
which is independent of multipole structure of the gravitational source. In a 
particular case when the electromagnetic radiation pulse is generated by the 
wave source during a finite time interval, the structure of the wave tail at the 
time after the direct pulse has passed the gravitational source is in the first 
approximation independent of the higher multipole moments of the source of 
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gravitation, including the angular momentum. These results are then applied 
to a compact binary system. It follows that under certain conditions the tail 
energy can be a noticeable fraction of the primary pulse energy, namely, it 
is shown that for a particular model the energy carried away by the tail can 
amount to 10% of the energy of the low-frequency modes of the direct pulse. 
The present results indicate that the wave tails should be carefully considered 
in energy calculations of such systems and that the delay effect of the wave 
tails may be of great importance for their observational detection. 
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Typeset using REVT^ 



I. INTRODUCTION 



The problem of radiation propagation on a curved spacetime has been studied since the 
birth of general relativity, both as a matter of principle and as a basis of observational 
predictions. A field satisfying a hyperbolic differential equation (wave equation) propagates 
not only on the characteristic surfaces (light-cones), but also inside the light-cone in the form 
of wave tails which violate the Huygens' principle [Q,^. Physically speaking, the wave tails 
arise because the radiation is backscattered by the spacetime curvature. In certain cases 
backscattering can influence observations as it weakens and disperses sharp initial pulses. 
For instance, the electromagnetic (or gravitational) radiation from a pulsed source in the 
vicinity of a massive body reaching an asymptotic observer is received as two distinct pulses: 
one arriving along the direct route and the other, the tail, effectively being scattered off by 
the central body 

Recently the wave tails have become to be recognized as factors in the planned obser- 
vational detection of the gravitational waves by forthcoming laser interferometric detectors 
||6HrT|. It has also been shown that the wave tails play an important role in the generation 



of gravitational waves by the orbital inspiral of a compact binary system |T^-|n| . The close 



relationship between the generation of the wave tails and gravitational focusing has been 
demonstrated in Ref. [113 . 



The process of wave propagation on a curved spacetime being quite complicated, usually 
the wave equation is solved in the weak-field and slow-motion limit by the method of succes- 
sive approximations, using multipole expansion in one way or another (see e.g. [p!^-[T9|). A 
general solution of the electromagnetic wave equation in a curved space was constructed by 



DeWitt and Brehme |^0[ in terms of the Green's function, using the Hadamard procedure 
The general solution demonstrates scattering of the electromagnetic radiation by the 
spacetime curvature. An alternative technique for investigating electromagnetic radiation in 
the Schwarzschild and Kerr metrics was worked out following the Regge- Wheeler approach 
to metric perturbations pl|-p3|. This technique relies on an expansion into generalized 
spherical harmonics and is especially useful when radiation emanates from a given multipole 
moment of the source. Herein, for a point charge an infinite sum emerges. The benefit, 
however, is that the obtained solution is also valid for the strong-field region. 

We have recently developed a new method P3-E^ for calculating the exact solutions 



of scalar and tensor wave equations whose source terms are arbitrary-order multipoles on 
a curved spacetime. The developed method is based on the higher-order fundamental so- 
lutions (Green's functions) for wave equations which are defined in our paper |^ as the 
distributions that satisfy wave equations with the corresponding order covariant derivatives 
of the Dirac delta function as the source terms. Provided that the classical Green's func- 
tion and the multipole expansion of the source term are given, there is no need for a small 
expansion parameter within the framework of our formalism and in certain cases we can 
even find the exact multipole solutions for a strong field. Our approach also enables us to 
find an approximate solution provided that an approximate form of the classical Green's 
function is known which is the case for most spacetimes. Moreover, it proves to be more 
advantageous to apply our algorithm p5| , |26| , instead of the traditional approach of succes- 



sive approximations, as the amount of computations involved is considerably reduced and, 
as will be demonstrated below, some features unrevealed by the successive approximation 
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methods are brought forward. It is worthwhile to point out that, as distinct from most of 
papers dedicated to the topics of wave tails in which the source of the gravitational field is 
regarded as a point mass, within the framework of our approach [0,^ the extension of the 
source of gravitation is finite. The last circumstance enables us to avoid in computations 
the additional non-physical singularity, the regularizing of which may bring on difficulties in 
interpreting the results. 

On the basis of our method of higher-order Green's functions we have developed a 
new approach to the electromagnetic radiation on a curved spacetime. In a recent communi- 
cation [^] we presented some initial results obtained within the framework of this approach. 
Namely, we considered a pulsed source of electromagnetic radiation in arbitrary bounded 
motion in a weak gravitational field and concluded that generally the received radiation 
tail arrives after a time delay which represents geometrical backscattering by the central 
gravitational source. This delay effect of the wave tails may be of great importance by 
their observational detection. Further, by applying this approach to a compact astrophysi- 
cal binary system we demonstrated that under certain conditions the tail energy can be a 
noticeable fraction of the direct pulse energy. The underlying formulas involved are herein 
published for the first time. 

The wider aim of the present paper is to provide a comprehensive view of our new ap- 
proach to the electromagnetic radiation, more fully describe the already published results 
and expand upon them. Specifically, among the rest we will prove the following, (i) Van- 
ishing of the Ricci curvature tensor is the necessary and sufficient condition for the validity 
of the Huygens' principle in the first-order approximation, as it is already known, (ii) If the 
direct pulse of electromagnetic radiation has passed the gravitational source, then in the 
first approximation the structure of electromagnetic wave tails is independent of the higher 
multipole moments of the gravitational source, including the angular momentum, (iii) In 
an arbitrary weak gravitational field it is valid that in the first approximation the nonlo- 
cal radiative electromagnetic tail term at infinity acquires a universal form, viz. Eq. (^S[), 
which is independent of the multipole order. (The limitations of applicability as well as the 
relationships of these results with the earlier ones will be discussed in due course below.) 

Finally, we think that our method deserves to be presented in detail as it may prove to 
be useful also for the theoretical investigations of detection of gravitational waves performed 
within the framework of the LISA mission. 

The remainder of the paper is organized as follows. Sec. II gives a review of the the- 
ory of classical and higher-order Green's functions (fundamental solutions) for vector wave 
equations, as well as the recurrent formulas for calculating the Green's functions proceeding 
from the Hadamard coefficients. In Sec. Ill we determine the multipole moments of electro- 
magnetic field source term with respect to a given worldline, and also present an algorithm 
for calculating the exact multipole solutions of the wave equations. In Sec. IV the main 
results are obtained. We consider the tail term of the retarded Green's function expanded to 
first order in the gravitational potential, and give the first-order tail term for the multipole 
solutions of electromagnetic field (IV A, B) . In Sec. IV D we turn to nonlocal radiative wave- 
propagation correction in the far wave zone. We estimate the magnitude of electromagnetic 
radiation energy and find that, compared to the energy of the direct pulse, the value of tail 
energy can be considerable and may thus have astrophysical significance. Limitations to the 
applicability of the derived conclusions are also considered. Sec. V contains brief concluding 
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remarks. Appendix explains our notation and displays the relevant definitions. 



II. METHOD OF HIGHER-ORDER GREEN'S FUNCTIONS 

To investigate the electromagnetic wave tails within the framework of general relativity, 
we first consider on a pseudo-Riemannian 4-space M a vector wave equation 

L u =f, (1) 

which in local coordinates can be written in the following coordinate invariant form (for our 
notation see Appendix) 

Lu, := g'^'VaV^u, - R^ua = fc, (r) 

where the contravariant components of the metric tensor g"''^ are assumed to be of dif- 
ferentiability class C°°, and Rac are the covariant components of the Ricci tensor. The 
inhomogeneous term f of Eq. (|l]) is in general a distribution, i.e. f G V^{Q). In order to be 
able to complete the construction of u, we restrict the solutions to a causal domain Q <0 M 
(see point 5 in Appendix and Refs. PJ29[|). 

The notation and basic definitions used in this paper will be in detail presented in 
Appendix. Here we touch on only the spacetime subdomains frequently resorted to. C^{y) 
is the future light cone, i.e. the set of all points x E fl that can be reached along future- 
directed null geodesies from y; C~{y) is the past light cone, defined similarly by past-directed 
geodesies. The sets D^iij) denote the respective interiors of the future and past light cones 
C±(|/), whereas J^{y) := D^{y) U C^{y). 



A. Classical (zeroth-order) Green's function 

The classical (i.e. zeroth-order) Green's function, or fundamental solution as mathemati- 
cians would say, G(x, y) of the wave equation (|I|) satisfies 

LG(x,t/) = g(x,?/)(5(x,?/), (2) 

where 5{x,y) is the Dirac delta distribution, with ((5(x, y), 0(x)) := 0(?/) for all G C^iVL) 
and g(a;, ?/) is a transport bivector (for its defining equations see Appendix). 

As in the case of fiat spacetime, there are two particularly important Green's functions 
of the wave equation (|l|): the retarded Green's function G"'"(a;, y) and the advanced Green's 
function G~{x,y). It has been demonstrated that they have the following form p|,|29[|: 



G^{x,y) = ^ [U(x,i/)4°^(a(x,y)) + V(x,i/)e±(a(x,y))] , (3) 

where the bitensors U G C°°(^7 x Q) and V G C°°{Q x Q) are the Hadamard coefficients of 
the classical Green's functions (|) of the vector wave equation (0). 

The bitensor U in Eq. (|^) is determined by the following transport equation and nor- 
malization condition: 
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where 



y)VaUl{x, y) + M{x, y)Ui{x, y) = 0, 

Ul{x,x) = 5l (4) 



M(x,y):=^V'^V,a(x,y)-4. (5) 



It is well known that the bivector U can be written as [29,20 



U := g(x,?/)A;(x,y), (6) 

where g(x, and k{x,y) are respectively the transport bivector and scalarized van Vleck 
determinant, defined in Appendix. 

The bivector V G C°°{Q x Q) in Eq. (|^), called the tail term, is determined by the 
characteristic Cauchy problem. It corresponds to the 'logarithmic term' of the Hadamard 
construction and is inherently connected with the concept and validity of Huygens' principle 
In the regions D^{y) the vector field V satisfies the homogeneous wave equation 

LV(x,y)=0, (7) 

which is completed by the characteristic initial conditions 

a^"(a;, y) V,n*(x, y) + (M(x, y) + 2)K^(x, y) = -hu^x, y), (8) 

for Vx G C^{y). 

Sometimes, instead of solving the differential equations (|^) and (|^) straightforwardly, in 
order to find the tail term V, it is preferable to use an exact integral equation. We proceed 
from the fact that Friedlander has derived the corresponding exact integral equation for the 
tail term of the Green's function of the scalar wave equation (Eq. (5.4.19) in Ref. [p9[|), 
and suggested a procedure for obtaining the tensor field Green's functions from the scalar 
Green's function. Thus, generalizing the above-mentioned Friedlander's equation, we find 
for the vector case when x G D~^{y), the tail term V satisfies the following exact integral 
equation: 



V:{x,y) + - V,^{x,z)LU;{z,y)fi^^,,y^{z) + - m{x,z)LU;{z,yM^^ (9) 

Zn J'S{y) ZtT JS{y) 

Here: the 3-surface S(?/) and the 2-surface S{y) are defined by E(?/) := C^{y) H J~{x) and 
Siy) := C^iy) n C~(x), respectively; the operator L acts at z; the 3-form iicj(z,y) and the 
2-form uj are defined in Appendix by Eqs. (|63|) and (|65|). 
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FIGURES 




FIG. 1. The geometric regions occurring in Eq. (9): an illustration on a 3-dimensional inter- 
section of the Minkowski space with the plane = const. The darker region on the future light 
cone C^iy) is the hypersurface whose 2-dimensional boundary S is represented by an ellipse 
which is the intersection of the light cones C~^{y) and C~{x). 

We will use Eq. in Sec. IV for obtaining an expression for electromagnetic wave tails 
in the first-order approximation. 

Only the retarded Green's functions G"*" will be discussed, as the corresponding results 
for the advanced Green's functions G~ can be obtained by reversing the time orientation 
on the domain f2. 



B. Higher-order Green's functions for vector wave equations 



Let f2 be a causal domain. Then the tensor differential operator L in (|lD has a unique 
retarded /ith-order Green's function G"'' on such that 



x,y) 



:-irv 



supp G^ c J^{y)- 

If test vector fields $ G V^{VL) are such that supp $ C then the retarded /ith-order 

Green's function G'*' is of the following form: 



G+ = — 
27r 



.v=0 



(10) 



where W(^, i/ = 0, . . . /x+ 1 are bitensor fields of rank 1 at x and of rank /i + 1 aX y recursively 
determined by 



W^i = W;:+2 := 0, W° := U, W? := V, 



(11) 
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where = 0, . . . , /i + 1. 

This statement was first proposed by us without proof in |27], its proof can be found in 
2^ . Eqs. ( ]Tl| ) provide a simple recurrent algorithm for calculating the /ith-order Green's 



functions for the vector wave equation ([^) proceeding from the Hadamard coefficients U and 
V of classical Green's function (^. The formulas (|TO|-|lT]) constitute the main mathematical 
tool which allows construction of our exact multipole solutions of vector wave equations. 



III. EXACT SOLUTIONS 
A. Multipole expansion of the electromagnetic field source 

On the basis of Dixon's ideas PD[ we will now construct a multipole expansion of the 



electromagnetic field source f as follows. We first choose a unique timelike worldline ^ 
lying inside the worldtube F C f2 of the source of the electromagnetic field that represents 
its dynamical properties. Such a curve can be given as a C°° embedding t ^ y(t) G ^ 
of an open interval I G TZ into fi, where TZ is the real line. We set v^{t) = this 
vector is assumed to be timelike and future-directed, and it is convenient to normalize the 
parametrization so that t is the proper time, which means that f*f « = 1. For y{t) E ^ we 
define := {x\v''{t)a.i{x,y{t)) = 0}, the spacelike hypersurface consisting of all geodesies 
through y(t) orthogonal to v(t). We suppose that there exists a 3-form /x^ such that fi{x) = 
drjA^s on supp f; here ri{x) is defined by v\ri{x))a-i{x,y{ri{x))) = 0, y{ri) G ^. For the sake 
of simplicity, we assume also that supp f is compact in the domain Q. Regarding the source 
function f as a regular distribution with compact support, we can write 



ff, 



dt < f(x)$(x) > IJ.T,{x) 
Jl JT,{t) 



where < f$ > denotes the scalar product. Let M^(i(:) be tensor fields of corresponding 
ranks /x + 1 > 1 at y(t) G ^, with supp M^(t) compact, and let A*" > be a integer. We 
consider the line distributions 

fN = E(-i)'^M^'^^^Wv^(,)^/^if (x,y(t))^?j^ G s'\n) (12) 

which assigns to any $ G C^{Q) the number 

N 

(f^., := E / Mi^''^Ht)^r,miymt- (13) 

In the particular case of = the formula (0) means that 

{gt{x,y{t))M^,{t)5{x,0,Mx)) := / M^(t)0.(y(t))rft. 



Let us now choose a test function so that supp $ D supp f . If for all test functions such 
that supp $ D supp f , the following equation is valid 
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(14) 



where $Af is determined by the Taylor expansion of the vectorfunction $ 

11=0 ^ f^- 

then the hne distribution is called the Nth-order multipole expansion of the source func- 
tion f . 



If is the A^th-order multipole expansion of f, then it follows form (|T^) and ([T^) that 
one can choose 



Ml^'^'it) = ^ j^^^^gi{x,y{t))f^{x)a'^^{x,y{t))...a'^'^{x,y{t))f^^{x). (15) 

The tensor field M^(t) determined by the expression (|15D is called the 2'^ -pole moment 
of the field source f with respect to ^. One should mention that Eqs. ([13|) and ([T^ ) do 



not uniquely determine the structure of the multipole moments in Eq. (0). For 
example, let us define a set of new multipole moments M^(t) by relations M^^^^^ (t) = 

Mf^^'^^it) + D^J'^^^^^v'^ + j-^Djj''^^ {t) , fi < N , D^i^^^ = , where the symbol | denotes 
the absolute derivative along the worldline ^. Now, integrating Eq. ([13D by parts, it follows 
that Tat remains unchanged for arbitrary tensor fields D^(t) with compact support. 

The multipole moments defined by Eq. (|15D are obviously symmetric in the first 
indices, i.e. M^i---*'^-? = M^*!---*^')-^. They also satisfy the orthogonality condition M^^ v,^ = 
0. One should say, though, that due to charge conservation, Va/" = 0, the multipole 
moments of different orders are interrelated through differential constraints. For that 
reason their direct usage in solving particular physical problems is inconvenient. However, 
Dixon demonstrated in Ref. that by means of the multipole moments a new set 



of reduced multipole moments , > 1 can be constructed, which together with the 
total charge Q will completely describe the four-current f. If Q = const, then the charge 
is automatically conserved, and thus the multipole moments of different orders , A* > 1 
are independent. Dixon's reduced multipole moments have the following symmetry and 
orthogonality properties: 

Qil^^'' = Q^Xi'"'^^'^ for /i>0, 

Q^:^^-'^^'"^ = 0, Q;^;r\, = o for /i>i. (16) 



Following Dixon's work [|30|, we will in this paper approximate the source of electromag- 
netic waves f by its iVth-order multipole expansion (O) replacing, however, the multipole 
moments M^(t) with the reduced ones M^(t): 

^ + 1 , r- - \ J 
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B. Exact solutions of vector wave equations with a multipole source term 



In one of our earlier papers we demonstrated how to obtain a solution by way of 
describing 2'^-pole radiation of a vector field in terms of the higher-order Green's functions, 
i.e. how to calculate the retarded solution u+ of the vector wave equation 

Lu^ = Pi 

with a multipole source term 

Pl := {-iYMl^^^^{t)VAUgZ)i^Mmi^^^^^^ (18) 

Here M^(t) is a tensor field of order /x+l on the worldline ^ of the source of (electromagnetic) 
radiation and the coordinates of the points of the line ^ are denoted by y{t), where the 
parameter t is the proper time along ^. Let us suppose that there is a to ^ such that 
= for t < Iq. If is determined by relations (0), it follows from the symmetry 
properties of Dixon's moment s (p!6D that VaP^ = 0, i.e. 

(v„p;:,0) = -(p;:,v.0) = o. (19) 

On the domain J'''(0\{0) with J^{i) := U?/e5 J^i.y)i ^^e solution u+ can be represented 
as a regular distribution (function). To find a form of the solution suitable for applications 
it is reasonable to define a new, retarded time coordinate r(x) as follows 

a(x,i/(r(x))) = 0, y{T)eC-{x). (20) 

Evidently, the future light cone C^{y{t)) is determined by the equation t{x) = t. We denote 
the corresponding Leray form by Pt{x), i.e. drAprix) = fi{x) on ^"""(OVIO' (^^e Ref. [^). 
On the surface C'^ {y{t))\{y{t)} we have 

da = -a.i{x, yit)y{t)\r(^)=tdT. (21) 

In what follows we shall use the notation 

^{x,t) := -a,i{x,y{t)y{t). (22) 

Due to the theorem 5 proved in [^, the unique retarded solution u+ of Eq. (pISD can be 
written in the following form: 



[X) 



1 ^ 

—y 



^XM-M/(^)^(t)H^,^;,(^)(x,y(t)) 



ip{x, t) dt J '^{x-, t) 

1 r{x) 



r(x) 



^2^17 ^r,m(^Mt))Ml^^^\t)dt, Vx e JHOm, (23) 

where the quantities W^, u = 0, . . . , fi are bitensor fields of rank 1 at x and of rank fi + 1 at y 
recursively determined by Eqs. ([TT|), and the bitensor field V is the tail term of the Green's 
function of Eq. (|l]). Identifying the tensor field in Eqs. (|l^) and (^31) with the tensor 
field (|13), defined via Dixon's multipole moments, it is easy to see that due to Eq. (|19|), 
and the uniqueness of the retarded solution of the wave equation, the solution satisfies 
the gauge condition VaW^" = 0. Consequently, the quantity m^" can be interpreted as the 



retarded potential of electromagnetic field. Thus Eq. (23) with the recurrence system (|TT]) 
enables one to find with admirable ease the exact multipole solution of arbitrary order for 
the electromagnetic wave equation (|18D by means of the world function a, transport bitensor 
K and tail term V of the classical Green's function 
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IV. ELECTROMAGNETIC WAVE TAILS IN CASE OF A WEAK 

GRAVITATIONAL FIELD 



Application of formulas (p3D presupposes knowledge of the tail term V of the classical 
Green's function and the world function a, the exact forms of which, however, have been thus 
far calculated for some particular cases only, such as for example the Bianchi-type 1 metric 



3T| , de Sitter metric and a class of Robertson- Walker metrics [^]. Nevertheless, for 
most cases approximate forms of V and a are known and can be used within the framework 
of our formalism in obtaining approximate solutions. 

We assume that the gravitational field is weak and therefore use a perturbational ap- 
proach. The metric tensor and other geometrical quantities are supposed to depend on a 
small parameter e which determines the order of deviation from fiat spacetime. Such expan- 
sion has been used for instance for the electromagnetic field in Ref. [0] and for gravitational 
waves in Ref. W^- 



A. First-order tail term for the Green's function 

We assume that the metric tensor gab is expanded up to the first order in e, so that the 
zeroth-order term is the metric tensor of the flat spacetime Qab^ i-e. 



9ab = gab + Slab + 0{e'^), 


gab ^ gab _ ^^ab ^ q^^2^ (24) 


with 7"* = g""^ g^'^'^cd- The world function a{x,y), the transport tensor g{x,y) and the 


scalarized van Vleck determinant k{x, y) can then also be expanded in the parameter e. We 
have 



a = a + ea + 0(e 
1 ^ 



2\ 



g = S + eg + 0{e^), 
01 

k = k + ek + 0{e^), (25) 
01 

where ex, g and = 1 are the corresponding quantities of fiat spacetime. The expansions 



of the Ricci tensor and the tail term V of the classical Green's function begin with terms 
which are small of the first order, that is 

Rab = £Rab + 0{e'^), 
1 

y = e\ + 0{e^). (26) 

Due to the tensor character of the expressions obtained, they can be used in all the coordinate 
systems, where the metric preserves the form (0). Up to now we have where necessary 
meticulously inserted the underscripts and 1 denoting the zeroth-order terms (i.e. the fiat 
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spacetime quantities) and the first-order terms, respectively. The Ricci tensor Rat and the 
tail term V being small of the first order, we will omit these underscripts where confusion 
is excluded, and write the results in what follows to the first-order approximation. 

Next we will find an explicit expression for the leading member of the expansion of the 
tail term V proceeding from the integral equation (^). Since V is small of the first order, 
it follows from Eq. that the quantity LUl{x, y) is also small of the first order. Thus the 
first integral in Eq. (^) is a second order small quantity, whereas the first factor Up{x,z) 
under the second integral in Eq. (|^) may be considered as a zeroth order quantity. Hence, 
taking into account Eq. (^) with k = I, the vector integral equation reduces to 

V:{x,y) = -^ j gl{x,z)LU;{z,y)uj{z) + 0{e'). (27) 

Z-K Js(y) 

Likewise, the first factor under the last integral is of zeroth order, and if we would use 
Minkowskian coordinates in the fiat spacetime we could write g'^ = 6^ + 0{e). Nevertheless, 
in order to be able to use general curvilinear coordinates in the fiat spacetime without 
any need for modifications in the integral in Eq. (p7|), we prefer to keep the more general 
expression gP{x,z). The same will be done in similar cases below. 

To calculate the factor LU^ under the last integral we need the linear terms in the 
expansions of ga{x,y) and k{x,y) which can obtained as follows. Expanding the defining 
equation a'"''Vagi{x,y) = of the transport bivector, we have 

ga^{x,y)= ga. + U\'a{y,z{X))gl{y,z{X))[^,^p,r){z{X))-^^^^^^^^^^ (28) 

^ "'0 '- 

Here, as well as in the next integral, integration is performed along the geodesic z{X) connect- 
ing the points x and y, where A is an affine parameter on the geodesic, such that z{l) = x, 
z{0) = y. 

Next substituting Eq. (||) with Eq. (|^) into Eq. (|) yields for the scalarized van Vleck 
determinant (0) the following first-order expression 

k{x, y) = l + la^'^ix, y)a-\x, y) [' A(l - A)(7^(x, ziX))glix, ziX))R,M>^))dX + Oie'). (29) 







Finally, after quite simple but voluminous calculations we obtain from Eqs. (p7| - p9D the 
seeked expression for the first-order tail term 

V:{x,y) = ^{9yb [cr-\x,y)P\x,y)]+a-\x,y)r,{x,y)}, (30) 



where 



P\x,y):= I gl{x,z)G^\z)di:,{z), (31) 
F^{x, y) := f gy'i URID.^ + RD^, + 2R„a{x, y)) uj{z). (32) 
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Here G^" : = 
Ricci scalar, 
defined by 



n 



pq 



\g^''R are the contravariant components of Einstein's tensor, R is the 
= a.^[p{z,x)a-^qT^{z,y) and dLqi^z) denotes the 3-surface element on E(?/) 



dz' A dz'' A dz' 



(33) 



where e^g^s are the components of the discriminant tensor with eoi23 = 1- 

It can be shown that V^{x, y) remains bounded when x — x G C^{y). From Eqs. pO| - |32D 
we see that vanishes for every x eVL and y EVtii Rpq{z) = 0, G Q. It is the necessary 
and sufficient condition for validity of Huygens's principle in the first-order approximation 
1 31,11 



Among the astrophysical applications of great interest is the case in which the source 
of gravitation is spatially isolated, i.e. supp Rab C T, where T is the world tube of the 
gravitational field source. In this case it is reasonable to divide the spacetime domain J~^{y) 
into the following three subdomains: 

D{y) :={x|a:G J+(y),S(y)nf = 0}, 
D*{y) := {x\x G C+iy) H f C S{y) H f = 0}, 



D := {x\x G J+{y), S{y) n f ^ 0}. 



(34) 



It is evident from Eqs. (pO|-p^) that V'{x,y) = if Vx G D{y). We remind that our 
results are valid for an arbitrary weak gravitational field. An analogous conclusion on 
the background of a weak Schwarzschild field was drawn in Refs. [p5|-|37[|. As in the first 
approximation P\ = 0, x G D*{y), the expression for the tail term V in the domain D*{y) 
becomes considerably simpler, namely 



V{x,y) 



1 

4^' 



{x,y)P'Vt 



(35) 



Thus, in the subdomains (|3^) the character of the tail term is qualitatively different, (i) 
No-tail region. In the subdomain D the observer at x does not see any wavetail, V(x, y) = 0, 
from the source at y. (ii) Simple-tail region. In the subdomain D* there appears a wavetail 
of the simple form (|35|) . (iii) General-tail region. In the subdomain D one can examine a 
wavetail of the general structure (PO) (see Fig. 2). 
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FIG. 2. A diagram in Minkowski 2-spacetime: the regions where the character of the tail 
term is quahtatively different, y is the source of electromagnetic waves and T is the world-tube of 
gravitational source. The observer in the region D(y) does not see any wavetail, the observer in 
D*{y) can see a wavetail of simple structure (37) whereas in the region D(y) one can see a wavetail 
of the general form (30). 

Eq. ( pSf ) leads to the conclusion that if the metric tensor satisfies the Einstein field 
equations, the tail term of the Green's function in the first approximation in the spacetime 
domain D*{y) is completely determined by the four- momentum of the source of gravitation. 
The higher multipole moments of the source of gravitation, including the angular momentum, 
do not influence in this approximation the structure of Y{x,y) (in the domain D*{y)). A 
corresponding result for a scalar field has been presented in Ref. p8|, and for a vector field 



On the basis of Eq. (p3D , the tail term of the retarded 2'^-pole solution of Eq. ([ISD 
can be written as 



Next we will study the tail term of a multipole solution in a weak gravitational field 
more closely. In a weak gravitational field the bivector field V(x,?/) is determined by Eq. 
(^). In what follows we presume that the worldline of the source of electromagnetic waves 
remains outside the source of gravitational field, whereas ?/(to) and y(ti) are points on the 
worldline ^ of the wave source, to and ti being, respectively, the proper time values when 
the source begins and terminates the emission. If the worldline ^ of the wave source lies 
outside the world tube T of the source of gravitation, i.e. supp C, fl supp T = 0, and if 
supp M^(t) C (to, ^i), then in the structure of the tail term V° there appear features similar 
to the case of the tail term of the fundamental solution V discussed at the end of the last 
subsection. 

To begin with, we divide the spacetime domain ^^(l/(io))\{0 with y{t) G C, into three 
sub domains: 



in Ref. H. 



B. Tail term for multipole solution 




(36) 



E := {x\x e J+(2/(to))\{0, s(y(to)) n r = 0} 



E* := {x\x e J+(2/(to))\{0, c+iyih)) n r c s(y(ti)), 

5(i/(ti))nf = 0}, 



E := {x\x G J 



iy{to))\{0,x^E,x^E*}. 



(37) 
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FIG. 3. A diagram in Minkowski 2-spacetime: the regions E, E* and E where the tail term of a 
wave pulse with finite duration behaves differently. The past light-cone C~ (x) originating from the 
point X is represented by the dotted lines. The bold-faced part of the future light-cone C~^{y(to)) 
corresponds to the hypersurface T,{y{to)) with x = x, and the boundary of the hyperplane S(y(to)), 
i.e., S{y(to)) is seen as the two points s. T is the worldtube of the source of gravitational field. ^ 
is the worldline of the wave source which radiates during a finite proper time interval [to,ii]. A is 
the worldline of the observer. On the observer's worldline A these intervals are indicated in bold 
where the observer can in principle see the direct pulse (the interval [a, b]) and the wave tail (the 
semi-interval [c, oo)). During the interval (6, c) on ^ there occurs a blackout between the direct 
pulse and the wave tail. 

In the subdomains E, E* and E the character of the tail term of a wave pulse with finite 
duration is qualitatively different, (i) No-tail region, li x & E, then also x e D{y{t)) for 
every i e (to, t(x)). Hence it is valid that — 0, and consequently 

V; = 0, Vx e E. (38) 

As E' 7^ 0, there exist points in space, where the (fore) front of the wave is not simultane- 
ously accompanied by the wave tail, the latter appearing after some time delay. To describe 
this effect in more detail, let us define the quantity To{x, z) as a solution of the following set 
of equations: 

cr{x,z) =0, z e C~{x) n f , 
a{z, y{To{x, z))) = 0, y{To) G C-{z). (39) 

For further treatment are important the maximal and minimal values of the time interval 
t{x) — tq{x,z) (measured in the proper time of the wave source) for the fixed spacetime 
point X, viz. 

Ai(x) := t{x) — max tq{x, z), 

A.2{x) :— t{x) — min To{x, z). (40) 

The quantities Ai(x) and /S.2{x) can be given a simple lucid interpretation by taking 
into account that an instantaneous wave pulse emitted by the wave source at max ro(a;, z) 
travelling at the speed of light reaches the spacetime point x exactly as if it would first travel 
to the source of gravitation, reflect from there and further travel to the spacetime point x. 
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In the case of min to{x,z) the reflection takes place precisely when the above-mentioned 
instantaneous wave pulse has passed the source of gravitation. 
Obviously, in case t{x) < to + ^i(^); then x E E and = 0. 

(ii) Simple-tail region. If t{x) > ti + A2(x), then x E E*, and the structure of the tail 
term is determined only by the four-momentum of the source of gravitational field: 



-[g]ix,yit))Ml^'^^^it)P'V, 



a{x,y{t)) 



dt. 



(41) 



Consequently, the higher multipole moments of the source of gravitation, among them 
the angular momentum, not influencing the structure of the wave tail. 

(iii) General-tail region. To the domain E corresponds the interval (to+ Ai{x) , ti-\- A2{x)) 
of the retarded time t(x), and for these instants it is valid 



min {T{x)—Ai{x),ti) 



= — 



V^^,^^){x,y{t))Ml^^)^{t)dt 



(42) 



to 



where V is determined by Eqs. ( |30| - p2D . Thus, in comparison with the direct pulse, the tail 
of the wave appears to the observer after the time delay Ai(a;). 

The particular case Ai(a;) > ti —t^ is also of interest. In it there occurs a time lapse of 
duration A(x) between the end of the principal pulse {t{x) = ti) and the appearance of the 
wave tail, with 



A{x) := Ai(a;)-(ti-to 



(43) 



Hence, instead of a single pulse, the observer would see two clearly separable pulses: the 
principal one and the wave tail. 



C. Schematic representation of the geometry of wave tail generation 

The purpose of this subsection is to provide a schematic plane spacetime representation 
of the geometry of wave tail generation. This subsection is meant as an illustration, the 
rest of this paper does not rely on it. For the sake of simplicity we make here the following 
simplyfing assumptions which are not used in our calculations: (i) the gravitational source 
is spherical and static, and (ii) the gravitational and wave sources, and the observer do not 
move with respect to each other. 
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FIG. 4. A schematic representation of the geometry of wave tail generation on a spacehke 
plane x^x^. Here y is the wave source, x is the observer, the darker shadowed circular disk is the 
gravitational source, and K is the focusing region. The remaining components in the figure are 
explained in Sec. IV C. 

In order to draw Fig. 4 we have suppressed the time dimension and one space dimension. 
Thus all the images in Fig.4 are obtained by projecting the Minkowski spacetime onto a 
hypesurface = const, and then cutting the 3-space with x^x^. The plane of the figure 
is determined by locations of the wave source y, of the observer x and of the center (not 
pointed out on the figure) of the gravitational source. The surface S{y), which spreads with 
time, is the boundary of the ellipsoid of revolution with foci at the locations of the 

observer x and of the wave source y. The ellipses 5*1 and 5*2 are the intersections of the plane 
of the figure with the surface S{y) that, respectively, correspond to the instants of time at 
which a delta-like wave pulse emitted by the wave source reaches and passes the source of 
gravitation. If the surface S{y) has not yet reached the gravitational source, then V = 
and there is no wave tail. If the surface S{y) has passed the gravitational source, the source 
will forever remain inside the region of integration while P{x,y) /Sir = const will be 

the 4- momentum of the gravitational source and F^{x, y) = 0. The shadowed area K, which 
includes the source of gravitational field, corresponds to the region where the tail wave field 
is predominantly generated by gravitational focusing which deforms the direct wave fronts 
||T5| . Thus we have the following picture. The wave source y emits an instantaneous wave 
pulse. The direct pulse propagates along the direct route yx to the observer x, after which 
there occurs a blackout before the arrival of the first tail contribution. The parts of the wave 
front which travel (along the routes h, I2, etc.) to the points of the region K, scatter off 
(reflect) from there and then propagate (along II, I2, etc.) to x appearing to the observer as 
the tail wave. 

To elucidate the introduced maximal and minimal time intervals, Ai and A2, let us again 
consider to the example presented in Fig.4. At the time y^ = to the wave source emits a 
delta-like pulse. The evolution of the surface S{y) is characterized by the time-varying semi- 
major a and semi-minor b axes which depend on the observer time x° as a = [x^ — to)/2 and 
b = (o"(x, ?/))^/^/2. The ellipses Si, S2 and 6*3 are the intersections of the plane of the figure 
with the surface S{y) at times of observation x^ = to + \x — y\ + Ai, X2 = to + \x — y\ + A2 
and Xg > X2, respectively. Here \x — y\ is the spatial separation between the wave source 
and the observer. The observer detects the direct pulse at the time Xq = to + \x — y\, the 
wave tail begins to appear at x^ = to + ^1 + ^1 = a^o + Ai, and beginning from the time 
X2 = to + I2 + I2 = + A2 the structure of the wave tail is determined by Eq. ( pO]) with 

= 0. The time interval Ai = li + II — \x — y\ is evidently equal to the difference of the 
propagation times of two pulses from the wave source to the observer: one arriving along 
the direct route and the other first travelling to the source of gravitation, reflecting from 
there and then reaching the observer. 

D. Radiative part of the electromagnetic wave tails 

Let us examine the case in which the source of gravitation is spatially isolated and the 
distance of the source of waves from the source of gravitation is bounded. The tail term 
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being a first-order small quantity, we can regard the spacetime as fiat when discussing it 
and use Minkowskian coordinates whose origin lies inside the world tube of the source of 
gravitation. We denote the distance of the observer to the origin of the coordinates by r, 
and will also below use the following notation: if two functions rA{x) and rB{x) are equal 
in the limit r oo, r = const, i.e. lim {rA{x)) = lim {rB{x)), then we will write 



r — 'OO 
T— const 



r — i-oo 
T— const 



A{x) = B{x). The position vectors of the points x and y in three-space are denoted by x 
and y, respectively. 

Taking now into account that for an arbitrary finite r are valid the asymptotic relations 
hm cr(x,y(t)) ~ r, lim a-Ax^yit)) ~ r and lim aik(x,y(t)) ~ 0(1), we obtain from Eqs. 

r— >o o r— >oo \ / / r— >oo ' 

-|3^) the relationship 



V;',U^){x,y{t)) = {-ira.,,{x,y{t))...a.,^{x,y{t)) 



d 



^ip{x^ t) dt ^ 



Vrix,yit)). 



(44) 



For brevity, in what follows we will not write explicitly the arguments of the functions 
depending both on the point x as well as on the parameter t. 

The radiative part of the wave tail which dominates at infinity can be written as 



rad 



t{x) 



27r 



cr-i, ... a 



to 



^ dt 



V^dt. 



(45) 



Due to Eqs. 
written as 



ja ( |TT| ) the zeroth-order radiative part of the 2''-pole solution can be 

1 d\'^ „ a-i, ... 0"-, 



Zl{xMr)) = 



= ^9]ix,yiT)) 



1 



2tt 

dV 



ip dt 



9 



-M^(^)^'(t) 



t=T{x) 



^^{x, t) dr ^ 



a.,,ix,yiT))...a.,^{x,yiT)) ^^^j^^y. ^ 
V'(a;,r) ^ 



(46) 



On the basis of the last two relations, after integrating by parts the right hand side of Eq. 
(^) we obtain the following expression for j.^^^: 



dt 



^fi rad 



t{x) 



to 



ijglVtZlixMmt 



-Ml}^'^^{t) 



(47) 



t=T{x) 



where Zj(x,?/(t)) := Zj(a;, y(r(a;)))|^(^.)=t. 

We must point out that Eq. (^) contains /i + 1 terms, out of which the last ones 
are actually instantaneous. The first term in Eq. (^T]) is truly nonlocal, whereas the factor 
ipglVf" assigns most of the weight to the source's recent past. 

The truly nonlocal radiative wave-propagation correction 



rad ■ 



t{x) 



ijglVtZl{x,y{t))dt 



to 



(48) 
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takes an universal form which is independent of the multipole order. An analogous result 
in a weak Schwarzschild field in a slow-motion approximation has been earlier found in Ref. 
p!0| (cf. also [0). Our formula ( ^HD generalizes this result in two ways. First, Eq. ( ^HD 
is valid in case of an arbitrary weak gravitational field in the corresponding asymptotically 
flat spacetime. Second, there are no formal restrictions to the bounded motion of the wave 
source, e.g. the wave source can move with a relativistic speed. Due to the universal form 
of the tail term, Eq. (^) has a simple, but from the point view of observations essential, 
generalization. On the assumption that the local part A"-{x,t) of a wave observable at 
infinity (or its zeroth-order radiative part) can be approximated with sufficient accuracy 
by the superposition of a finite number of multipole waves, it follows from Eq. ( ^HD that 
the corresponding nonlocal radiative wave-propagation correction E"-{x, r) has the following 
form: 

E'^ix, t) = r^^^ iljglA\x, t)Vtdt, (49) 

•I to 

where ip is defined by Eq. (p^ . 

For a source radiating in the pulsed mode the last formula enables us to comparatively 
simply (on the basis of the parameters of a recorded direct pulse and a presumable model of 
an astrophysical object) predict the physical parameters of the radiation tail and evaluate 
the possibilities of observational detection of the tail. 

It should be emphasized that applicability of Eq. ( P?[ ) does not depend on whether the 
source of radiation actually emits in the pulsed mode or pulses of radiation are present at 
the location of an observer due to the kinematics of the radiating system (e.g. the rotation 
of the directed radiation cone of a pulsar). 

If A^(a;, r) does not vanish merely during a finite time interval tf) < t < ti and the 
worldline of the source of radiation lies outside the world-tube of the source of gravitation, 
then the conclusions of Sec. IV B are also valid for the radiative part E°- of the tail term. 
Thus, for example, if r(x) < to + ^1(2^)5 then E°-{x, r) = 0; if r(x) > ti + A2(x), then 

E^{x, r) = r>(a;, ^x, y{t))P'V, ( ^ ) dt 



1 



At: J to a{x,y{t)) dt 

A matter of particular interest is the situation in which Ai{x) > {ti — to), as in this 
case there is a blackout during the time interval A{x) = Ai(x) — (ti — tg) between the 
end of the direct pulse and the appearance of the wave tail. The last fact considerably 
simplifies the observation methods for distinguishing the profile of the direct pulse from 
the general relativistic radiation tail originating from compact astrophysical binary systems. 
The relative intensity of the direct pulse and the wave tail, and the time delay (blackout) 
between them can yield essential additional information, independently of other methods, 
about the physical characteristics of a binary (the distance between the compact objects, 
the mass of the source of gravitation, the orientation of the plane of the orbit with respect 
to the observer, etc.). 
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From Eq. (^) it follows that the wave tail effect of the astrophysical systems radiating 
in the pulsed mode is predominatly caused by the low-frequency modes (ti; ^ 27r/(ti — to)) 
of the direct pulse. (Here the angular frequency u of the radiation is measured in the proper 
time of the wave source.) The high-frequency waves reflecting from the spacetime curvature 
interfere and prevalently attenuate each other in the expression of the tail. Also important, 
from the point of view of observational detection, is the fact ensuing from Eq. (|49D that the 
intensity and the time delay of the wave tail of radiation from a pulsed source moving on a 
circular orbit depend substantially on the mutual positions of the observer, the wave source 
and the source of gravitation. By comparing the profiles of the pulses emitted in different 
points of an orbit, one can distinguish in principle the contributions of the direct pulse and 
the tail even if there is no considerable time delay between the tail and the principal pulse. 
This circumstance is significant because in the case of most astrophysically realistic models 
the physical conditions necessary for the occurrence of a blackout in between the direct pulse 
and the wave tail will considerably decrease the intensity of the tail. 

Let us now turn to the magnitude of electromagnetic radiative energy arriving after the 
direct pulse has passed, t(x) > ti + A2{x). The power ^ radiated into the solid angle dil 
can be calculated as follows 

dP 

^^^^ ^dn= -^^:o(^'^)^a,o(x,r). (51) 



For further treatment it is reasonable to present Eq. (|50|) in a somewhat different form. In 
the proper reference frame of the source of gravitational field we have 

El ^ -2M ' l^°°'f';'f'"-'r'°'d«, (52) 

Juito) (n(r) - u)2 

where u{t) := y^(t) — n ■ y{t), n := x/\x\ and M is the mass of the source of gravitational 
field. Below we will use notations T := u{ti) — u{tQ), which is the duration of the direct 
pulse in the observer time, and A2 := u(ti + A2{x)) — u(ti). 

To get an idea of the magnitudes involved we construct an artificial example in which 

A'^^{x,Tix))\r=tiu) = ^r(t;,0)e(M-M(to))e(M(ti)-M). (53) 

Here G is the Heaviside distribution, v, (j) are the polar angles determining the observer's 
position and is an arbitrary vector function. From Eqs. (|5Tl-p^) we obtain an estimate of 
the ratio of the intensity of the wave tail J(ti + /S.2{x)) to the intensity of the direct pulse 
Jo (the ratio of the densities of the energy fluxes at the observer's location), namely 



J(ti + A2(x)) 




(54) 



We see that if A2(x) is sufficiently small (A2(x) ~ 2M), the intensity of the tail can be 
comparable with that of the direct pulse. Evidently, in this case we cannot confine ourselves 
to the first approximation but must instead additionally take into consideration the higher 
approximation(s) . 

To illustrate the above estimate, let us consider a wave source rotating in a circular orbit 
(of radius tq) around a spherical source (of radius r^) of gravitational field. The origin of the 
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spatial coordinates is taken to coincide with the center of the source of gravitation. Under 
the conditions n ■ y(ti) < and tq ^ d + Vg, where d := rQ^l — {ft ■ y(ti)/ro)^, it follows 
from Eqs. (M) and {mj that 



(d + Ts 

2rn 



+ 0(M). 



(55) 



Hence, if tq > r^/4M, then in case the wave pulse is emitted in the region of the geometric 
shadow of the source of gravitation or in its vicinity it is valid that A2(x) ~ 2M, and 
according to the estimate (^4]) the intensity of the tail can be of the same magnitude as the 
intensity of the primary pulse. 

For the model under discussion it is also easy to find the ratio of the energy S transferred 
by the tail term (beginning from the time ti + A2(a;)) to the energy of the direct term Sq, 
namely 



T 



where 



Fip) =P + 



2p^ln.^ + ' 



p+ 1 p 

It is interesting to note that the function F[p) has a maximum at p ~ 0.638. 




P. =0.638 1 



FIG. 5. The function F{p) vs its argument p. 

Therefore, for the given model the amount of energy transferred by the tail after the 
time ti + A2{x) is maximal if the duration of the direct pulse is T ^ 1.57A2(a;), and can in 
this case make up nearly 10% of the energy of the direct pulse: £ ~ 0.119£^o(2M)^/A2(a;). 

As the quantity JF does not explicitly depend on the mass of the gravitational source, then 
its value for any particular case (as 0.119 for the present example) is primarily determined 
by the profile of the direct pulse and the spatial configuration of the system consisting 
of the wave source, gravitational source and the observer. Somewhat unexpected is the 
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outcome that the magnitude of the factor (2M/A2)^, which characterizes the influence of 
the gravitational source, can be of the order 1, even if the potential of the gravitational source 



is low everywhere, i.e. 2M -C r^. This conclusion can be understood on the basis of Ref. |T5 
where it is shown that the tail field is predominantly generated by the direct field in those 
regions where gravitational focusing has deformed the geometry of the direct wave fronts (see 
the shadowed region in Fig. 4). The deformation is characterized by the focusing function 
a{z,y). If the points z and y lie on a geodesic line which does not cross the gravitational 
source, then a{z,y) = 0. If the distance rg of the wave source from the gravitational source 
is much larger than the extension of the latter, then for the rays originating in the wave 
source and passing through the gravitational source it is valid a ~ Mro/r^. In the case 
of our example the condition 2M/A2 ~ 1 corresponds to a ~ 1. Let us note that the 
last effect is not revealed by the traditional methods based on the expansion in terms of 
spherical functions, as in this case it is natural to choose the effective location of the wave 
source within the source of gravitation. Thus, A2 ~ and a ~ M/r^, that is, the intensity 
of the tail is proportional to the square of the potential of the gravitational field on the 
surface of the source of gravitation. 

Similar estimates for the energy of scalar and electomagnetic dipole wave tails in a strong 
Schwarzschild field have been found respectively in Refs. p9|,E0|. 



E. Constraints on applicability of the elaborated algorithm 
Relying on Ref. ITH], we will now briefly analyse limitations to applicability of the above 



results. Firstly, note that the conclusions of this paper were derived under the assumption 
that on the physically interesting spacetime domain the world function a is single- valued, 
i.e. the geodesies originating from the point y do not cross. If this condition is not fulfilled, 
then the developed algorithms are not directly applicable. However, by virtue of the su- 
perposition principle, valid because of linearity of the wave equation, these algorithms can 
be applied even when the world function is a multiple-valued function. Then the correct 
classical Green's function is the sum over all distinct elementary classical Green's functions 
corresponding to all distinct geodesies between the points x and y (see also pO|). 

Crossing of geodesies would be caused by gravitational focusing; and at the crossing 
point the exact factor k{x,y) in U (see ( |2^)) would diverge. Thus, the criterion for no 
crossing is that k{x,y) be finite along C^{y). Let us now consider our first-order expression 
(p9|). This expression for k{x,y) can never diverge if the gravitational source is bounded. 



However, if the focusing function a{x,y) = k{x,y) — 1 approaches unity, then the second- 
and higher-order effects will come into play, producing a divergence. Thus, the constraint 
a{x,y) < 1 for all x and y is necessary, on the one hand, for the first-order analysis to be 
valid; one the other hand, it simultaneously avoids crossing of geodesies. 

For a wave-source in a circular orbit (the observer lying on the plane of the orbit), Eq. 
(PPP gives max a ~ {2M/rs){2rQ/rs), where M is the mass of the source of gravitation, 
is its linear size and tq is the radius of the orbit. In this case the constraint max a < 1 is 
significant: it says that in order to avoid too much ray focusing, the wave-source and the 
gravitational source must not be too far apart 

^ < ^. (56) 
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Applying Eqs. p9|), (|50D and (0) to a particular model system, one should not overlook 
the fact that, although E is a first-order small quantity and therefore, when calculating 
the integrals involved the spacetime may be considered as a flat background; generally the 
retarded time r(x) must be regarded as the retarded time on the curved spacetime because 
the integrands may be sensitive to the wave phase. 

The quantities A2(x) and Ai(a;) := u(to + — u{to) offering their own inherent 

physical interest must generally also be calculated at least up to the accuracy of the first- 
order corrections. 

Let us now consider in more detail the time intervals A2(x) and Ai(x) for the binary 
system described at the end of Sec. IV C. We denote the origin of the coordinates at the 
center of the source of gravitation by z and consider only the case in which the points x, y 
and z are aligned on a straight line. Evidently there are two possible cases: (i) the point z 
is placed between the points x and y, (ii) the point y is placed between the points x and z. 
Equations and (^D|) yield the following conclusions. In the case (i) the time intervals 
Ai(x) and A2{x) have maximal values: 

max Ai(x) = 2(ro - r^) + 



max A2(x) = 2(ro + r^) + X2{x). 
In the case (ii) the time intervals Ai{x) and A2(x) have minimal values: 

min Ai(x) = 0, min A2(x) = + Xsfx). 

2ro 

Here the quantities Xa ~ 2M, a = 1,2, 3, denote the corrections caused by the Shapiro 
effect within the inner region of the binary system which are of the same order of magnitude 
as the Shapiro time delay between the points y and z. 

It is obvious that in Eqs. (^), (|50|) and ( ^2l) the corrections caused by the inner Shapiro 
effect may be neglected only in the low-frequency limit 

The necessary condition (^) is also sufficient, if the inequality (|56|) is simultaneously satis- 
fied. 



V. SUMMARY AND DISCUSSION 



i) On the basis of our earlier works [P7| , p4| , p5| , we developed simple recurrent formulas, 



i.e. Eq. (|TT|), for calculating exact multipole solutions, i.e. Eq. (^), of the electromagnetic 
wave equation on the background of a curved spacetime which proceed from the classical 
Green's function (i.e. the corresponding Hadamard coefficients). 

(ii) The next main new contribution of this paper is formula (|3^) with ([30| ) for the first- 
order tail term of the electromagnetic multipole wave in case of an arbitrary weak background 
gravitational field. From these formulas will follow a number of physically interesting effects. 
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Firstly, one can infer from Eqs. (pOD the earlier known (see p^ , |33[] ) necessary and suffi- 
cient condition for the validity of the Huygens' principle (in the sense of Hadamard) in the 
first approximation, namely Rat = 0. 

Secondly, if the gravitational field source is spatially isolated, then beginning from a 
certain instant of time the structure of the tail accompanied by a pulse of electromagnetic 
radiation is completely determined by the four-momentum of the source of gravitation, and 
the higher multipole moments, including the angular momentum, do not at all infiuence the 
structure of the tail in the first approximation. For the scalar wave equation and for the 
Green's function of the electromagnetic wave equation essentially similar conclusions have 



been presented in Refs. |Q and p3| . 

Thirdly, in Sec. IV B we discussed a further consequence of Eq. (^), namely, the time 
delay effect of the tail with respect to the direct pulse. In the case of compact binary astro- 
physical objects the time interval (blackout) in between the direct pulse and the tail caused 
by delaying of the tail may be observable and can with the relative intensity of the tail give, 
independently of the other methods, essential information about the characteristics of the 
physical system (the distance between compact objects, the mass of the source of gravitation, 
the orientation of the plane of the orbit with respect to the observer, etc.). Earlier on the 
basis of the structure of the Green's function of the wave equation in a weak Schwarschild 
field the time delay effect of the tail and its possible astrophysical implications have been 
indicated in the papers and [Q. Let us mention that some unpublished calculations by 
the present authors demonstrate the occurrence of a similar delay effect also for linearized 
gravitational waves in the second approximation (cf. [^). By the common approach in which 
the wave equation is solved by separating the angular variables the delay effect, as a rule, 
remains unrevealed. This is due to the circumstance that ensuing from the symmetry of 
gravitational field (Schwarzschild, Kerr etc.), the origin of the spherical coordinates whose 
worldline is simultaneously the worldline of the multipole radiation source lies inside the 
source of gravitation where the Ricci tensor does not vanish. Rat 7^ 0. 

(iii) It is valid in the first approximation that the nonlocal radiative wave-propagation 
correction at infinity takes a universal form (^8|) which is independent of the multipole order. 
A corresponding result in a weak Schwarzschild field in a slow-motion approximation has 
been published in |]10|- Our Eq. (|4^) generalizes their result in two aspects. First, Eq. (^Sf ) 
is valid in case of an arbitrary weak gravitational field in the corresponding asymptotically 
fiat spacetime. Second, there are no formal restrictions to the spatially bounded motion of 
the wave source, e.g. the motion of the wave source can be relativistic. 

(iv) The intensity of the electromagnetic wave pulse emitted by the wave source within 
a compact binary in the vicinity of the geometric shadow of the source of gravitation can 
be of the same magnitude as the intensity of the direct pulse, and the energy carried away 
by the tail can amount to 10% of the energy of the low-frequency modes of the direct pulse. 
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APPENDIX: NOTATION AND BASIC FORMULAS 



1. We adopt the metric signature (+, — , — , — ) with the Latin indices running and sum- 
ming from to 3. To abbreviate the notation of repeated tensor indices, we introduce 
multi-indices, e.g. A{fi) := (oi . . . a^), /(/i) := (ii . . . i^), -P(/i) := {pi ■ ■ -Pfi), etc. We apply 
extensively the technique of two-point tensors (or bitensors) and follow closely the standard 
notation and index conventions To distinguish the tensor indices which refer to the 
point X from those referring to the points y and z we use indices a,b,. . . A{fj,) , B (fj,) , . . . 
at X, indices /(/i), . . . a.t y and p,q, . . . , ... at Thus F/j-g^ is a 
contravariant vector at x and third-order covariant tensor at y. The index convention is also 
used for ordinary tensor fields, in which case it distinguishes the value of the component v°' 
of a vector field v at x from its value at y, and G'/(^) denotes the covariant components 
of the tensor field G of rank /i at the point y. 

2. Symmetrization and antisymmetrization of the tensor indices is denoted respectively 
by the parentheses and brackets. For example := {f"^ + t''")/2 and tt'^''! := (t"^ - t''")/2. 

3. Ordinary differentiation is denoted by a comma (,). Covariant differentiation 
with respect to the Levi-Civita (metric) connection is denoted by V and semicolons, e.g. 
^x-f{x,y) = f;aix,y), g^^W^bfi^x.y) = /'"(x,?/) and V/(2)/(a;, ?/) = \/y^2'^y^lf{x,y) = 
f-^i^i^{x,y). Absolute differentiation along a line z{\) is denoted by an overdot, i.e. 

= w]jZ^ . 

4. We use a system of units in which the speed of light and the gravitational constant 
are equal to 1. 

5. The class of C°° tensor fields of rank m is denoted by £'^[VL) and the subspace of 
£^■{0) consisting of fields with compact support by T''"(f2). A tensor- valued distribution 
T{x,y) G P'™'(r2) of rank m dX x E Vt and of rank at ?/ G is a continuous linear 
map 'D'^{Vt) —>■ S^{Q). If {uj,tt) is a coordinate chart such that x,y & uj and $ G T>"^{uj), 
then each component of {TfQ^\x,y), (f)A{m){x)) is a (scalar- valued) tensor distribution (for 
a detailed discussion see |^). A set fio C is called past-compact if J^{x) nf2o is compact 
(or empty) for all x E Q. We denote the class of distributions in T>'{Q) with past-compact 
supports by V^{Q). 

6. The world function biscalar a{x,y) G C°°{yL x Vt) is defined by the equations 

g'^'Vacyix, y)Vi,a{x, y) = Aa{x, y), (58) 
a{x,y) = a{y,x), limcr = 0. (59) 

The world function cr{x,y) is numerically equal to the square of the geodesic arc length 
between the points x and y, and is positive for timelike intervals and negative for spacelike 
intervals. 

7. The propagator of geodetic parallel displacement (also called the transport bitensor) 
is defined as a bitensor field g(x, y) of rank m at both x and y, which satisfies, in local 
coordinates, the following differential equations and initial conditions: 

^^"V,4':;))(x,y) = 0, (60) 
9V)('^,^) = K---^lZ- (61) 

8. The surface distributions 6^\a{x, y)), jj, = 0, 1, 2, . . . are defined by 
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{6i^\a{x,y)),<P{x)) ■= ^^^J'Pl' / , (62) 

for all G C^{Q), if /i = 0, and for all G if > 1, where /icr(a;) is a Leray 

form for which 

da{x,y) A fia{x) = fi{x), (63) 

fi{x) being an invariant volume element and Cf{y) = {x\x G D^{y),a{x,y) = 6,e > 0}. In 
general, the integral (p^ ) can be evaluated by means of a partition of unity subordinated to 
a covering of Q by open coordinate neighbourhoods in each of which 



/i(x) = \Jg{x)dx^ A dx^ A c/x^ A dx^, (64) 

where (^(x) := \det{gab{x))\ is the determinant of the metric tensor. The relevant properties 
of the surface distributions are given in p9|j41[] . 

9. The 2-form uj on the 2-surface S{ii) := C^{y) H C^{x) in the integral (|]) is defined by 

da{z,x) /\ da{z,y) /\uj{z) = jj,{z), (65) 

where /i(z) is an invariant volume element on M. 

10. The Heaviside distributions are defined by 



(e±(a(x,?/)),0(x)) := / 0(x)/i(x), (66) 

where J^{y) '■= D^{y) U C'^ly). The relevant properties of the Heaviside distributions can 
be found in Ref. |2|gl[. 

11. The scalarized van Vleck determinant is a biscalar which is defined by the equation 

,/ N \det a,ai{x,y)\^^'^ , . 

M9{x)9{y)\^'* 

12. A nth-order Green's function of the tensor wave operator L with respect to [y, Q) is 
a tensor- valued distribution G{x,y) G of rank 1 at x and of rank fi + 1 aX y which 
satisfies the equation 

LGl^,)ix,y) = i-irVBi49t{x,y)9^l;:^ix,y)5ix,y)] , (68) 
where is the tensor propagator of geodetic parallel displacement. 
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